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We discuss a semiclassical calculation of low energy charge transport in one-dimensional (Id) 
insulators with a focus on Mott insulators, whose charge degrees of freedom are gapped due to 
the combination of short range interactions and a periodic lattice potential. Combining RG and 
instanton methods, we calculate the nonlinear ac conductivity and interpret the result in terms 
of multi-photon absorption. We compare the result of the semiclassical calculation for interacting 
systems to a perturbative, fully quantum mechanical calculation of multi-photon absorption in a Id 
band insulator and find good agreement when the number of simultaneously absorbed photons is 
large. 



I. INTRODUCTION 

In Id electron systems, interactions have a very pro- 
nounced effect and give rise to a variety of unusual 
phenomenal. Interacting Id systems with a lattice po- 
tential display insulating behavior at commensurate fill- 
ing. Examples include charge density waves (CDWs), 
and fermions on a lattice with local Hubbard interac- 
tions. Due to interaction effects, the spin degrees of free- 
dom are gapped and decoupled from the charge degrees 
of freedom, such that they do not influence the dynamics 
of charge transport. In these systems, the charge degrees 
of freedom are described by the quantum sine-Gordon 
model. 

The low frequency response of a half filled system 
with periodic potential is characterized by an optical 
gap. Hence at zero temperature the linear ac conductiv- 
ity vanishes for frequencies smaller than the gap, apart 
from a series of discrete resonances which appear for suf- 
ficiently strong repulsive interaction^!. For frequencies 
larger than the gap energy the linear conductivity shows 
a power law dependence on frequency^. Similarly, the 
linear dc conductivity shows a power law dependence on 
temperature T at temperatures larger than the spectral 
gap. 

At zero temperature and frequency, charge transport 
is only possible by tunneling of charge carriers, which 
can be described by instanton formation. The nonlinear 
dc-conductivity is characterized by££ / ~ exp(—E /E), 
and for fixed value of the electric field strength is a mono- 
tonically increasing function of the ac frequency of the 
external field. In this work, we revisit the instanton cal- 
culation of the nonlinear ac conductivity by Maki and 
discuss in detail the physics of multi-photon absorption in 
the dynamic limit. In addition, we point out that in the 
static limit the efficiency of soliton-antisoliton pair an- 
nihilation and the dynamics of scattering between these 
quasiparticles decisively determines the value of the con- 
ductivity. 

To be specific, we consider the quantum sine-Gordon 
model. We first scale the system to its correlation length, 
where the influence of the potential is strong and a semi- 
classical instanton calculation becomes possible. We re- 



produce an earlier result^ for the nonlinear ac current 
within the Matsubara formalism, and carefully analyze 
the dependence of the instanton creation rate on the ra- 
tio of electric field energy to photon energy. In the dy- 
namic limit, the result can be interpreted in terms of 
multi-photon absorption. We compare the result of our 
semiclassical analysis to a fully quantum mechanical per- 
turbative expression for a Id band insulator. When ex- 
pressing the semiclassical result in terms of the energy 
gap, the solitonic correlation length and effective veloc- 
ity, it can be compared to the quantum mechanical ex- 
pression. In the limit where the photon energy is much 
smaller than the optical gap, there is good agreement 
between the two expressions. 



II. PERTURBATIVE CALCULATION FOR 
NONINTERACTING ELECTRONS 

As a reference point for the instanton calculation for 
interacting systems, we start by describing a fully quan- 
tum mechanical calculation of the conductivity of non- 
interacting Id electrons with a periodic potential at half 
filling. The periodic potential splits the band into a va- 
lence and a conduction band, such that the system be- 
comes a band insulator. When exposing the system to 
a monochromatic electric field E(t) = £<] cos(u;t), the 
only contribution to the ac conductivity comes from the 
excitation of electron-hole pairs. Clearly, the linear ac 
conductivity vanishes as long as the photon energy hw 
is smaller than the band gap -E gap - In the following, we 
sketch a derivation of the nonlinear ac conductivity. 

We follow the derivation described by Wherrett£. We 
consider a situation in which N photons are needed to ex- 
cite an electron from the valence to the conduction band, 
that is Nhiu > E gap . The coupling to the electromagnetic 
field has its origin in the A ■ p term in the Hamiltonian, 
its matrix elements are given by 
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In the last step it was assumed that the initial and final 
momentum are close to the Fermi momentum, and that 
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the ratio pF/m e g = v c g defines an effective velocity. We 
calculate the transition rate from Fermi's golden rule as 
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is 



given by p(e) 



The density of states 

(y / hk F V c ff/€)/(2V2TThv cS ). 

The transition amplitude Mvc^ between valence and 
conduction band involves N — 1 virtual intermediate 
states. For odd N, these states correspond to the con- 
secutive excitation and deexcitation between valence and 
conduction band. Assuming that the transition matrix 
elements have no important energy dependence, it can 
be calculated as 
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(2.3) 



In the last step, it was assumed that E g ss Nftw and Stir- 
ling's formula for the factorial function was used. Putting 
everything together, the transition rate per unit length 
is found to be 



For K > K c (u) the potential is RG irrelevant and de- 
cays under the RG flow, while for K < K c (u) the poten- 
tial is relevant and grows. For the periodic potential in 
the action Eq. (|3.ip . one finds K c (0) — 8/p 2 . We assume 
K < K c (u) and scale the system to a length £ = ae l , on 
which the potential is strong. After the scaling process, 
the parameters K , v, and u in Eq. (|3 . 1 [) are replaced by 
the effective, i.e. renormalized but not rescaled, param- 
eters if off, Vctt, and M cff . 

The compressibility k = is used as a generalized 
density of states for interacting systems, it is given by 
K c ff = Keff/vegir. Our calculations are valid for photon 
energies hu> and electric field energies eEo^K c s below the 

soliton energy E s = , 8 ■> ■> . In this RG calculation, we 

do not attempt to treat a possible nonlinear dependence 
of coupling parameters on the external electric field. The 
full inclusion of the external field in an equilibrium the- 
ory is not possible as it renders the ground state of the 
system unstable. The quantum sine-Gordon model has 
an infinite number of ground states connected by a shift 
of the phase field by 27r. Here, we concentrate on renor- 
malizing each of these ground states separately and take 
into account the coupling between different ground states 
due to the external electric field in the framework of an 
instanton approach. 
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From this absorption rate, the iV-photon contribution to 
the nonlinear conductivity can be obtained via 
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III. MODEL 



IV. INSTANTONS 

We consider a time dependent external field E(t), 
which upon analytical continuation itv e ff —> y turns into 
a field E(y). In imaginary time, the electric field has 
to obey the same periodic boundary condition E(y + 
0Ve«) = E(y) as other bosonic fields, e.g. the displace- 
ment field f(y). This boundary condition is respected by 
a discrete Fourier representation^ 



E(y)=Tj2E("n)e 
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(4.1) 



The charge degrees of freedom of interacting Id elec- 
trons subject to a periodic potential are described by the 
quantum sine-Gordon model 
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where we have rescaled time according to vt — > y, and 
(3 = X/ksT. The dissipative part of the action describes 
a weak coupling of the electron system to a dissipative 
bath, for example phonons. It is needed for energy relax- 
ation in soliton-antisoliton annihilation processes^. We 
assume it to be so small that it does not influence the RG 
equations for the other model parameters significantly. 
The smooth part of the density is given by —d x (p, and 
p = 1, 2 for CDWs and LLs, respectively. 



with Matsubara frequencies to n . A monochromatic exter- 
nal field is hence described by E{y) = E cos(ui n y/v e s). 
Analytic continuation from Matsubara frequencies to real 
frequencies is defined by ito n — > lj + irj. The original 
calculation^ did not make use of Matsubara frequencies 
and used a dependence on imaginary time 



E(y) = E cosh (uiy/v eS ) 



(4.2) 



which does not respect the periodicity requirement in 
imaginary time and grows exponentially in y. We define 
the ratio 
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which is proportional to the ratio of the field energy ac- 
quired on a length scale -Kcff£ an d the photon energy hw. 
The static dc limit corresponds to 7 — > 00, whereas the 
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(a) 




FIG. 1: Instantons for different ratios 7 = & ° E ff cff ^ of field 
strength to frequency, (a) for 721, the instanton is almost 
circular, and (b) for 7 = 0.2, the instanton is elongated in 
time direction. 



dynamic ac limit corresponds to 7 — > 0. It turns out that 
for 7 > 1 it does not matter whether one analytically con- 
tinues to real frequencies before solving the equation of 
motion for the instanton or in the very end of the calcula- 
tion. However, in the range 7 < 1, the two prescriptions 
lead to different results. Using an oscillatory time depen- 
dence as defined by a monochromatic E(u n ) in Eq. (|4.ip . 
the instanton solution becomes periodic in time and has 
no well defined beginning or end. On the other hand, 
using the time dependence Eq. (14. 2p in the equation of 
motion, the instanton is well defined and the final re- 
sult agrees quantitatively with the quantum mechanical 
calculation presented in the previous section. Although 
it is not very intuitive why the exponential time depen- 
dence Eq. (|4.2p should describe an electric field periodic 
in time, it leads to the correct result and we adopt this 
prescription in the following calculation. 

The wall width 1 / y'Mefi ~ £ of an instanton solu- 
tion to the action Eq. (|3.1|) is for weak external fields 
much smaller than the extension E s /eoEo of the instan- 
ton. The domain wall energy per unit length is given 
by E s /v e g, and in contrast to the disordered case^, 
the spontaneous formation of instantons due to quan- 
tum fluctuations needs not be considered as the domain 
wall energy is spatially constant in the present model. 
Hence, the instanton action can be expressed in terms of 
the domain wall position X(y) as 
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The spatial extension of the instanton is 

r 1 \ 2Vctf 4. 1 
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The limiting forms for the static limit 7 
dynamic limit 7 — > are 



(4.6) 
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In the static limit, the instanton is circular with an ex- 
tension inversely proportional to the electric field. In the 
dynamic limit, the scale for the instanton extension is set 
by the period of the ac electric field, and it is elongated 
in time direction by a factor ln(2/7). 

The creation rate of kink antikink pairs is proportional 
to the negative exponential of the instanton action di- 
vided by the space-time volume of an instanton. One 
finds 
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This result would be changed by a preexponential factor 
when including the contribution of fluctuations of the 
instanton shape, of zero modes, and of the Jacobi de- 
terminant originating in the transition from an integral 
over the field tp to an integral over the the domain wall 
positio n 11 ^ 2 . The dependence of the exponent on 7 is 
displayed in Fig. [2j For the static limit with frequencies 
ui <C eo-EoC; 7 is large and we find 
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According to the expansion Eq. (|4.10p . the first correc- 
dyE(y)[X + (y) — A_(y)] (4.4) tion of order w 2 reduces the static value of the exponent. 



The equation of motion can be integrated to give 

sinh (ojz/v e ff) 



Hence, for a fixed field amplitude Eq, the creation rate 
increases monotonically with frequency. 

In the dynamic limit of large frequencies, one finds 



X±(y) = T7 / dz—= 



Jl — 7 2 sinh 2 (uz/Ven) 



(4.5) 
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The domain of integration is bounded by the singulari- 
ties ±j/o of the integrand and the solution corresponds in- 
deed to an instanton with finite Euclidean action (Fig.[T]). 



This result gives rise to an instanton creation rate 
J / e Q E K cS £,p eir 
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Scaling parameter y 

FIG. 2: Dependence of the exponent of Pe on 7 = 
e ° g fta/ fcff Z ' ^' ie ex P onen t grows monotonically with 7 and 
saturates for 7 — > 00. 

This instanton creation rate can be interpreted in terms 
of multi-photon absorption: the optical gap in the system 
is 2E S , and hence the smallest integer larger than the 
ratio jr^j- determines the number of photons necessary to 
create one soliton-antisoliton pair. The ratio of electric 
field strength to photon energy in the brackets agrees 
exactly with the noninteracting result Eq. (|2.4|) if one 
identifies E g = 2E S , uses the definition of the soliton 
energy, and considers the fermionic case p = 2. 

The present semiclassical approximation can be ex- 
pected to work well for hu <gC E s when the photon num- 
ber is large and the electrical field can be treated un- 
quantized. However, it does not reproduce the threshold 
behavior expected whenever the photon energy is tuned 
through an integer fraction of the optical gap 2E S . The 
specific physics of a Mott insulator as compared to a band 
insulator is contained in the energy dependent form fac- 
tor, which for the band insulator is just the density of 
states. The energy dependence of the form factor is prob- 
ably contained in fluctuation corrections around the sad- 
dle point, which would be interesting to calculate. The 
main conclusion from the agreement between the semi- 
classical and the quantum mechanical result is that the 
dominant electric field dependence for the Mott insulator 
is similar to that of the band insulator when expressed 
in suitable effective quantities, which parameterize the 
interaction strength. 

V. NONLINEAR CONDUCTIVITY 

The quantum sine-Gordon model is integrable, and in 
the absence of additional interaction terms or a coupling 



to a dissipative bath its solitonic excitations have an in- 
finite life time. According to this logic, even an infinites- 
imally small density of solitons would give rise to an in- 
finite dc conductivity. This conclusion does not sound 
realistic, and in the following we will discuss the behavior 
expected from physical systems which are approximately 
described by this model. 

As both kinks and antikinks carry a topological quan- 
tum number, they cannot simply decay but only anni- 
hilate each other 13 . Pairwisc annihilation gives rise to a 
decay rate Pdccay ~ niq and hence n ~ yj P(Eq,u) in 
equilibrium. As kinks and antikinks are solitary waves, 
this annihilation is only possible in the presence of a dis- 
sipative bath. 

In the dynamic limit with multi-photon absorption, the 
kink and anti-kink created by the electric field stay close 
to each other and it seems plausible that they will an- 
nihilate quickly such that the nonlinear ac conductivity 
can be calculated according to Eq. (|2.5|) . 

The static limit is more subtle. Maki££ suggested to 
calculate the conductivity according to Eq. (|2.5p in the 
dc limit as well, which would imply that every kink and 
antikink that come close to each other annihilate with 
certainty. This assumption would only be justified if 
the coupling to the dissipative bath was strong. If the 
coupling to the dissipative bath was weak, there would 
be a finite density of kinks and antikinks, and the con- 
ductivity would depend on their mean free path. Under 
the assumption^ that the mean free path is of the order 
of the inter-particle distance, the product of equilibrium 
density and mean free path would be of order one, and 
the conductivity would not be exponentially suppressed 
anymore, a result quite different from the exponential 
suppression in the case of strong dissipation. 

In summary, we have discussed the nonlinear ac con- 
ductivity of interacting Id electron systems in a periodic 
potential. We found that the dominant physical process 
in the dynamical limit is multi-photon absorption, and 
that the result of an instanton calculation is in quantita- 
tive agreement with a fully quantum mechanical calcu- 
lation for noninteracting electrons in the limit of a large 
number of simultaneously absorbed photons. 
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